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Abstract 

A QCD phase transition in the early universe could have left inhomo- 
geneities in the baryon to photon ratio and in isospin that might have affected 
nucleosynthesis later on. At very high temperature QCD plasma can be de- 
scribed by perturbation theory because of asymptotic freedom, but a possible 
phase transition requires a nonperturbative approach like lattice gauge theory. 
Assuming that a first order transition did occur, a dynamical set of equations 
can be solved to evolve the universe through it and to quantify the scale of 
inhomogeneity. Unfortunately this scale appears to be too small by two orders 
of magnitude to affect nucleosynthesis. 



1 Introduction 

Who has not wondered what happens when matter gets squeezed to higher and 
higher densities or heated to higher and higher temperatures? At first molecules 
would dissociate into atoms, then atoms would ionize, then atomic nuclei would 
undergo a liquid-gas type phase transition to a vapor of protons and neutrons. At 
high baryon density there will be so many protons and neutrons per unit volume 
that their constituent quarks will no longer know to which nucleon they belong. 
This is expected to happen when nucleons are about a factor of two closer to each 
other than they are in atomic nuclei, that is, at a density of about 8-10 times that 
in a nucleus. (Normal nuclear density is 0.15 nucleons per cubic fermi - equivalent 
to 2 x 10 14 g/cm 3 .) At this point quarks will be free to move over distances much 
greater than the confinement distance of 1 fm. This type of matter is a cold quark 
gas. It is possible, but not yet known, that it exists in the cores of neutron stars. 

At high temperature there will be many particle-antiparticle pairs produced. 
The lightest hadrons, the pions, will be the first to be produced in abundance, 
quickly followed by all the hadrons listed in the Particle Data Tables PJ. Hagedorn 
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realized already in the 1960's that the level density of hadrons grows exponentially 
dN/dm ~ exp(m/To) where m is the hadron mass 0. The Hagedorn temperature 
To has the numerical value 160 MeV. A thermal ensemble of hadrons treated as 
noninteracting particles therefore has a limiting temperature equal to T because 
of this exponential growth. In a sense it is the end of the hadronic world; higher 
temperatures bring a quark-gluon plasma. Intuitively, this transition comes about 
because the hadron density is so great that the constituent quarks and gluons are 
confused about which hadron they ought to be confined to. 

The main interest in a cosmological hadron to quark-gluon plasma phase tran- 
sition arises from its potential to influence big bang nucleosynthesis, discussed in 
section 2. Whether or not QCD with its known set of parameters undergoes a first 
order transition or something smoother is still not settled. This issue is discussed in 
section 3 using perturbation theory and lattice gauge theory. Assuming that there 
is a first order phase transition one needs nucleation theory to understand how the 
transition proceeds, and this topic is discussed in section 4. The cosmological QCD 
phase transition is then analyzed in section 5. Conclusions, and paths for future 
work, are given in section 6. 

2 Inhomogeneous Big Bang Nucleosynthesis 

A cosmological first order phase transition at T ~ 160 — 180 MeV could leave 
spatial inhomogeneities in the baryon to entropy ratio and in the ratio of protons 
to neutrons. If these inhomogeneities survive to T ~ 0.1 — 1 MeV they could 
influence nucleosynthesis. This was first pointed out and analyzed by Witten ||, 
by Applegate, Hogan, and Scherer M, and by Alcock, Fuller, and Mathews H. 
In thermal and chemical equilibrium one might expect that the baryon density in 
the quark-gluon phase is higher than in the hadron phase. This is called the baryon 
density contrast. Assuming a critical temperature of 160 < T c < 180 MeV, Olive and 
I computed this baryon density contrast to be 1.5 to 2.5 when hadronic interactions 
were neglected and 5 to 7 when they were included ||. One would expect the 
last regions of space to undergo the phase conversion to contain more baryons per 
unit volume than the first regions to phase convert because of the lack of time for 
baryons to diffuse. After phase completion the neutrons will diffuse more rapidly 
than protons because they are electrically neutral and therefore do not Coulomb 
scatter on electrons. This leads to isospin inhomegeneities, at least temporarily. 

A detailed calculation of inhomogeneous nucleosynthesis with a comparison to 
observed abundances of the light elements was done by Kurki-Suonio, Matzner, 
Olive, and Schramm [0. They considered baryon density contrasts ranging from 1 
to 100 and a fraction of matter in the high density regions ranging from 1/64 to 1/4. 
The average separation of the high density regions I was left as a free parameter 
as was the average baryon to photon ratio of the universe. Differential diffusion of 
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Figure 1: Conservative upper limit to the baryon to photon ratio 77 from the 4 He 
abundance Y p < 0.248 and the deuterium abundance D/H> 1.5 x 10 -5 . The thick 
curves are for volume fractions covered by the high density regions of 1/(2^2) (solid), 
1/8 (dashed), and 1/(16 v^) (dot-dashed). The thin curves are for 1/64 (solid) and 
1/256 (dashed). From |TJ. 



protons and neutrons was accounted for and then a standard nucleosynthesis code 
was run. By fitting the observed abundances of 4 He, D, 3 He, and 7 Li they concluded 
that the baryon to photon ratio must lie between 2 x 10~ 10 and 7x 10~ 10 (or 2 x 10~ 10 
if certain constraints on 7 Li were relaxed). They also concluded that / < 150 m at 
the time of nucleosynthesis. At the completion of the QCD phase transition this 
upper limit would have been about 1 m. Early estimates of this distance were 
uncertain but promising || |9|. See also the reviews listed in flQ|| . A better estimate 
of this scale is the purpose of sections 4 and 5. 

Recently the inhomogeneous nucleosynthesis calculation was redone with tech- 
nical improvements and updated estimates of the cosmic abundances of the relevant 



light elements by Kainulainen, Kurki-Suonio, and Sihvola [11]]. Their results are 
shown in figure 1. The high density matter was distributed in spheres. The in- 
homogeneities are ineffective in influencing nucleosynthesis unless the high density 
regions were separated by more than about 150 m at T = 1 MeV. 
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3 QCD Equation of State 



The QCD equation of state can, in principle, be computed from the functional 
integral formula for the partition function |L2] . 



J [d$\[diP}[dA<$[dC a }[dC a ] exp ifdrj d 3 xC eff \ (1) 



Here tp is a quark field, is the gluon field, and C a is the ghost field. The La- 
grangian is 

C eff = £ - i + /j^/ + 9^C 8 + <?/^, W , (2) 

C = -\FrF^-^ f (i@-m f -g4 a G a )i>f, (3) 
Fr = d»A v a -&>A»-g3 ahc A»A v c . (4) 

The subscript / refers to quark flavor. The G a are the generators of the color group 
SU(N C ). The quark fields are Grassmann variables antiperiodic in imaginary time 
r = it, < t < [3 = 1/T, the gluon fields are periodic in this time, and the ghost 
fields are Grassmann variables periodic in imaginary time. (The complex ghost field 
may be thought of as subtracting 2 of the 4 degrees of freedom implied by A^ to yield 
the 2 physical transverse degrees of freedom of the massless gluon.) The Lagrangian 
is written in covariant gauge with p the gauge parameter. Physical observables are 
independent of p. 

There are two well-known ways to evaluate this functional integral. One is 
perturbation theory in the coupling constant g; the other is numerical estimation 
using Monte Carlo techniques on a discrete lattice. Before reviewing results obtained 
in these approaches, let us try to predict what the result might be. Suppose that all 
quarks were massless and that all chemical potentials were zero. Then there is only 
one parameter in the QCD Lagrangian, the coupling g which is dimensionless. If 
this coupling is really a constant then dimensional analysis implies that the pressure 
must be some function of g times the temperature to the fourth power: P = f(g)T . 
Then the entropy density is s = 4f(g)T 3 , and the equation of state may be written 
as e = 3P, the same as a free gas of massless particles. This is a trivial equation 
of state without a phase transition. In contrast the quantum version of the theory 
displays dimensional transmutation. The dimensionless coupling constant g gets 
replaced by a renormalization group running coupling. At finite temperature and 



at one-loop order []DJ 



2 -2 = 247r2 

9 ^ 9 (llN-2N f )\n(T/A QCD ) 1 ' 

for an SU(N C ) gauge group (N c = 3 for QCD) and Nf flavors of quark. The scale 
parameter Aqcd, which has dimensions of energy, enters as a result of the necessity 
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to regulate and renormalize, that is, to render correlation functions finite at some 
arbitrarily chosen energy or momentum point. As a result the equation of state with 
P = f(g(T))T can now become complicated and not predictable by dimensional 
analysis alone. This allows nontrivial thermodynamics and the possibility of a phase 
transition. This is on account of quantum mechanics, not classical interacting field 
theory. 

The partition function for QCD was calculated within perturbation theory to 
order g 3 in the late 1970's by myself and then to order g 5 in the mid 1990's by 
Arnold and Zhai [ 15|] , by Kastening and Zhai [pl| , and by Braaten and Nieto [17] . 



For Nf flavors of massless quarks and zero chemical potentials the pressure is 



P = & ^T' 
45 
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1 + — N f 1 + — N f — 

32 ' 4 V 12 V 7T 

3 / 2 f a (M\\ zl2 f n (M\\ 2 
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a, 



(6) 



Here a s = g 2 /4n, and _F 4 and F 5 are lengthy untransparent functions of the number 
of flavors and of the renormalization scale M. The latter should be chosen judiciously 
so as to minimize the contribution of higher order terms in the expansion. For sake 
of illustration consider Nf = 3, M = 2ttT, T = 500 MeV, and a s ~ 0.2. Then 

P = T 4 [1 - 0.18 + 0.30 + 0.06 - 0.37...] (7) 

36 

representing the contributions of a s to the power of 0, 1, 3/2, 2, and 5/2. At so low 
a temperature, only about three times the expected T c , this expansion is apparently 
not very quickly converging. Although very useful for very high temperature because 
asymptotic freedom forces the running coupling to zero logarithmically, perturbation 
theory cannot reliably predict the location or order of a possible phase transition. 

By its very nature perturbation theory using free particle basis states is incapable 
of describing a phase transition. An alternative is to discretize space and (imaginary) 
time and estimate the functional integral expression for Z numerically with Monte 
Carlo techniques JTBJ. One works with a spatial volume V = (N s a) 3 and a time 



dimension (3 = N T a, where a is the lattice spacing and N s and N T are the number 
of sites in each direction. The thermodynamic limit is a — > 0, iV s — > oo, iV r — > oo, 
subject to ft held fixed and V — ► oo. The difficulties associated in reaching for the 
thermodynamic limit include the enormous computer power and storage necessary 
for the huge number of degrees of freedom inherent in a quantum field theory, and 
the theoretical difficulty with describing fermions on a lattice. A typical lattice for 
pure gauge theory with no quarks is 64 3 x 8, and with quarks is 16 3 x 4. These are 
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Figure 2: The equation of state of pure SU(3) gauge theory with no quarks. From 
0- 

not very large lattices, especially in the latter case. Nevertheless computer power 
and storage is increasing with time, and practitioners are constantly improving the 
numerical algorithms. For example, improved lattice actions, where derivatives are 
better approximated via higher order finite differences and short distance behavior 
can be supplemented with perturbation theory on account of asymptotic freedom, 
have helped significantly. Here I have relied on a recent review by Karsch ||19|| . 

Pure quarkless SU(2) gauge theory is expected to undergo a second order phase 
transition reflecting color deconfinement at high temperature whereas SU(N C ) with 
N c > 2 is expected to undergo a first order phase transition. The SU (3) equation of 
state computed in lattice gauge theory is shown in figure 2. One physical quantity 
is required to set the overall energy scale. This is chosen to be the string tension a 
at zero temperature; its numerical value is taken to be the same as in our physical 
world. Then the critical temperature is about 270 MeV and the latent heat is about 
0.5 GeV/fm 3 . Below T c is a gas of (presumably) weakly interacting very massive 
glueballs. The lightest one is 1700 MeV. Above T c all 16 degrees of freedom of 
massless (interacting) gluons are liberated. These computational results seem to be 
stable with time and therefore well established and accepted. 
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Figure 3: T c as a function of the lightest pseudoscalar meson (pion) with both 
measured in units of the lightest vector meson mass (p or u) calculated by several 
groups with a variety of lattice methods. From ||19|| . 

A compendium of results from different groups for QCD with two degenerate 
flavors of quarks is shown in figure 3. In the limit that the quarks are massless a 
second order transition is expected, with no proper thermodynamic transition when 
they are massive. Even with massive quarks one sees a sharp rise in the energy 
density in a narrow window of temperature. The middle of this range can be used 
to define a crossover temperature still called, somewhat incorrectly, T c . In order to 
set the physical energy scale the mass of the lightest vector meson, the p-meson, is 
calculated and set equal to its real world value of 770 MeV. The mass of the lightest 
pseudoscalar meson, the 7r-meson, is calculated and compared to the real world value 
of 140 MeV. The pseudoscalar mass is a monotonically increasing function of the 
quark mass so it can be used to determine what quark mass ought to be used in 
the lattice calculation to represent as closely as possible the real world. The T c is 
plotted against the ratio of the pseudoscalar to vector meson masses, in units of 
the vector mass. Taking the limit of zero quark mass is notoriously difficult. The 
actual range of computed pion masses starts at 1.2 GeV and decreases to 300 MeV 
but no lower. There is a lot of scatter among the results from different groups using 
different methods but the trend is towards T c = 170 MeV in the zero quark mass 
limit. 
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Figure 4: Energy density in units of T 4 as a function of T for two light quarks (2 
flavor), two light quarks and one heavy quark (2+1 flavor), and three light quarks 
(3 flavor). From P[. 



QCD with three flavors of massless quarks is expected to undergo a first order 
phase transition. The question is: What is the order of the transition for the real 
world with up and down quark masses of order 5-7 MeV and a strange quark mass 
of order 110-150 MeV? Some results are shown in figure 4. The lattice is 16 3 x 4 in 
size. The quark masses are actually temperature-dependent (the usual situation for 
numerical reasons) with mii g ht = 0.4T and m^ eavy = T. Results for two light flavors, 
two light and one heavy flavor, and three light flavors are shown in this figure. The 
energy density e, in units of T 4 , is plotted against T/T c . In all three cases there is 
a big jump in the energy density over a narrow window of temperature. The 2+1 
flavor case is almost indistinguishable from the 3 flavor case in the critical region. 
The statistical error bars are very small. More significant are the systematic errors 
that are much more difficult to estimate. Due to the relatively small lattice size and 
still rather large light quark mass it is difficult to say whether or not the real world 
has a first order phase transition. Even if not, the large jump in energy density may 
make the real situation practically the same as if it did. 
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4 Relativistic Nucleation Theory 



The dynamics of first-order phase transitions has fascinated scientists at least since 
the time of Maxwell and Van der Waals. Much work on the classical theory of nucle- 
ation of gases and liquids was carried out in the early part of the 20th century. There 
were and still are many important applications, such as cloud and bubble chambers, 
freezing of liquids, and precipitation in the atmosphere. The goal of nucleation the- 
ory is to compute the probability that a bubble or droplet of one phase appears in 
a system initially in the other phase near the critical temperature. Homogeneous 
nucleation theory applies when the system is pure; inhomogeneous nucleation the- 
ory applies when impurities cause the formation of bubbles or droplets. For the 
application we have in mind, namely the early universe, it seems that homogeneous 
nucleation theory is appropriate. In the everyday world it is usually the opposite; 
dust or ions in the atmosphere are much more efficient in producing precipitation. 
Nucleation theory is applicable for first-order phase transitions when the matter 
is not dramatically supercooled or superheated. If substantial supercooling or su- 
perheating is present, or if the phase transition is second-order, then the relevant 
dynamics is spinodal decomposition. 

Suppose that a system is cooled below its critical temperature. Then there exists 
a critical sized droplet (or bubble, depending on whether the energy density in the 
lower temperature phase is greater or less than the higher temperature phase). If 
a droplet that forms because of statistical fluctuations is too small, its surface free 
energy is relatively large and the cost in total free energy is positive. The droplet 
will evaporate. If the droplet is large, its surface free energy is unimportant, and the 
droplet will accrete molecules and grow. A droplet of critical size is metastable, it 
is balanced between evaporation and accretion. The classical theory of Becker and 
Doring p0 |, which is nicely reviewed by McDonald plf , says that the probability 



per unit time per unit volume to nucleate the dense liquid phase from a dilute gas 
is given by 

I = <Q (^P) 1/2 n(l)e-^ , (8) 

where e(i*) is the formation energy of a critical sized droplet consisting of molecules, 
prime denotes differentiation with respect to the number of molecules i, n(l) is the 
density of single molecules and a{i*) is the accretion rate of single molecules on a 
critical droplet. Usually the accretion rate is taken to be 

a{Q = -n{l)vA-Krls , (9) 

which is the flux of particles (v is the mean speed of gas molecules) striking the 
surface of the critical droplet times a 'sticking fraction' s less than one. The first 
term in this formula is a dynamical factor influencing the growth rate, the second 
term characterizes fluctuations about the critical droplet, and the product of the 
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third and fourth terms gives the quasi-equilibrium number density of critical sized 
droplets. 



The modern theory of nucleation was pioneered by Langer p2[ . Langer's theory 
is based in a more fundamental way on the microscopic interactions of atoms and 
molecules. It can also be applied close to a critical point where in fact most of the 
current interest in the condensed matter community has been. Langer's modern 
theory of nucleation yields the following formula for the rate: 

/ = f fi e- AF * /T (10) 

Z7T 

where AF* is the change in the free energy of the system due to the formation of 
the critical droplet. Q is a statistical prefactor which measures the available phase 
space volume. 
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Here £ s is the correlation length in the gaseous phase, k is a dynamical prefactor 
which determines the exponential growth rate of critical droplets which are per- 
turbed from their quasi-equilibrium radius r*. 

K = -ln[r(t)-r«,] (12) 

The basic structure is the same as in the classical theory, but the prefactors are 
different. The dynamical prefactor has been calculated by Langer and Turski [^3[ p4j] 
and by Kawasaki |25| for a liquid- gas phase transition near the critical point, where 



the gas is not dilute, to be 

2XaT 

K ~¥nfl- (13) 

This involves the thermal conductivity A, the surface free energy a, the latent heat 
per molecule i and the density of molecules in the liquid phase rig. The interesting 
physics in this expression is the appearance of the thermal conductivity. In order 
for the droplet to grow beyond the critical size, latent heat must be conducted away 
from the surface into the gas. 

For a relativistic system of particles or quantum fields which has no net conserved 
charge, such as baryon number, the thermal conductivity vanishes. The reason is 
that there is no rest frame defined by the baryon density to refer to heat transport. 
Hence the Langer formula cannot be applied to such systems. The extension to such 
systems was made by Csernai and me p6[ . The change in free energy due to the 
appearance of a bubble of hadronic matter in quark-gluon plasma is 

Air 

AF = — r 3 [P q (T) - P h (T)] + 47rrV , (14) 
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where r is the radius. The critical sized bubble has radius 

2a 



P h (T) - P q (T) 

which leads to 



(15) 



AF* = — arl . (16) 

The expression for the statistical prefactor is unchanged, but the formula for the 
dynamical prefactor, or growth rate, is proportional to the shear rj q and bulk ( q 
viscosities in the quark-gluon plasma instead of the thermal conductivity. 

_ 4a(3( q + 4r] q ) 

3(Aw) 2 r3 {U) 

This is inversely proportional to the square of the enthalpy (w = e + P) difference 
between the two phases. The nucleation rate is 

4 / ay/i g(3( q + 4 Vq )r , AR/T 
I= n{3f) 3(Awye q ' (8) 

For numerical purposes we use an MIT bag model type equation of state with 
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P q = (45.5 + 14.25) — T 4 - B 

P h = (5.5 + 14.25) ^T 4 . (19) 

The 45.5 approximates the effective number of degrees of freedom arising from mass- 
less gluons and up and down quarks and a strange quark mass comparable to the 
temperature. The 5.5 approximates the hadronic equation of state near T c arising 
from a multitude of massive hadrons. The 14.25 arises from photons, neutrinos, 
electrons, and muons common to both phases. The bag constant B is chosen to give 
T c = 160 MeV. Furthermore Csernai and I [|27| estimated that o = 50 MeV/fm 3 , 
£g = 0.7 fm, and r] q = 18T 3 p8| . Generally for relativistic systems the bulk viscosity 
is small in comparison to the shear viscosity and so we neglect it. 



5 Cosmological QCD Phase Transition 

Given the nucleation rate one would like to know the (volume) fraction of space 
h(t) which has been converted from the quark-gluon plasma to the hadronic gas at 
the proper time t in the early universe. This requires kinetic equations that use the 
nucleation rate I as an input. Here we use a rate equation first proposed by Csernai 
and me [0. The nucleation rate I is the probability to form a bubble of critical 
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size per unit time per unit volume. If the system cools to T c at time t c then at some 
later time t the fraction of space which has been converted to the hadronic phase is 



h(t)= dt'I(T(t'))[l-h(t')}V(t\t). (20) 

V(t', t) is the volume of a hadronic bubble at time t which was nucleated at the 
earlier time t'; this takes into account bubble growth. The factor 1 — h(t') takes into 
account the fact that new bubbles can only be nucleated in the fraction of space 
not already occupied by the hadronic gas. This conservative approach neglects any 
spatial variation in the temperature. However, it is an improvement over the formula 
proposed by Guth and Weinberg [29] for first order cosmo logical phase transitions 



in that it allows for completion of the transition: the Guth and Weinberg formula 
is only valid during the early stage of the transition. 

Next we need a dynamical equation which couples the time evolution of the 
temperature to the fraction of space converted to the hadronic phase. The following 
analysis is very similar to the one used to study the electroweak phase transition in 
the early universe by Carrington and me |3(| . We use Einstein's equations as applied 
to the early universe, neglecting curvature. The evolution of the energy density is 

— = -— , (21) 
dR R ' v ; 

where R is the scale factor at time t. This assumes kinetic but not phase equilibrium, 
and is basically a statement of energy conservation. We express the energy density 
as 

e = he h (T) + [l-h}e q (T), (22) 

where €h and e q are the energy densities in the two phases at the temperature T. 
There is a similar equation for the enthalpy w. The time dependence of the scale 
factor is determined by the equation of motion 



1 dR 9,-wGe 

Rii = i—- (23) 

This expression can be used to relate the time to the scale factor using the normal- 
ization R(t c ) = 1. 

We also need to know how fast a bubble expands once it is created. This is 
a subtle issue since by definition a critical size bubble is metastable and will not 
grow without a perturbation. After applying a perturbation, a critical size bubble 
begins to grow. As the radius increases, the surface curvature decreases, and an 
asymptotic interfacial velocity is approached. The asymptotic radial growth velocity 
will be referred to as v(T). The expected qualitative behavior of v(T) is that the 
closer T is to T c the slower the bubbles grow. At T c there is no motivation for 
bubbles to grow at all since one phase is as good as the other. The bubble growth 
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Figure 5: Temperature as a function of scale factor. 



velocity was studied by Miller and Pantano [ 31 1 . Their hydrodynamical results may 
be parameterized by the simple formula 



vj — 3 ( 1 



rp . 3/2 

TJ 



(24) 



which indeed has the expected behavior. A simple illustrative model for bubble 
growth is then 

V{t\t) = f (r.(r(f)) + Jl dt"v{T{t"))f . (25) 

This expression can also be written in terms of R, R', R" instead of t, t', t" . 

We now have a complete set of coupled integro-differential equations which must 
be solved numerically. These equation take into account bubble nucleation and 
growth, energy conservation, and Einstein's equations. They make no assumption 
about entropy conservation. 

Figure 5 shows the temperature as a function of the scale factor. For practical 
purposes, nucleation begins near the bottom of the cooling line. Thereafter, nu- 
cleation and growth of bubbles releases latent heat, which causes the temperature 
to rise. The increasing temperature shuts off nucleation, and the phase transition 
continues due to the growth of already nucleated bubbles. The temperature can 
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Figure 6: Average bubble density as a function of scale factor. 

never quite reach T c ; if it did, bubble growth would cease and the transition would 
never complete. This is a result of the equations of motion and is not an imposition. 
Figure 6 shows the average bubble density 

n(R(t)) = f dt'I(T(t'))[l - h(t')\ (26) 
Jt c 

as a function of the scale factor. The bubble density rises rapidly just before R 
reaches 1.007 and reaches its asymptotic value just after 1.007. 

Figure 7 shows the nucleation rate as a function of scale factor. The rate has 
a very sharp maximum between 1.0070 and 1.0071. The turn on and turn off of 
the nucleation rate corresponds precisely with the fall and rise of the temperature 
shown in figure 5. 

Figure 8 shows the fraction of space h which has made the conversion to the 
hadronic phase. When h = 1 the transition is complete and the temperature will 
begin to fall again. This occurs when Rf « 1.4464, to be compared with the value 
one would obtain from an ideal Maxwell construction -RMaxwcii = (239/79) 1//3 = 
1.44630.... In fact the whole curve h(R) is very close to the ideal Maxwell construc- 
tion, apart from its delayed start apparent in the figure. The interested student will 
work out the Maxwell formula from the equations given here. 
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Figure 9: Average bubble radius as a function of scale factor. 

Figure 9 shows the average bubble radius as a function of scale factor. 

Air 

-^-f 3 n = h (27) 

It grows with time and with the scale factor, of course. At the end of the phase 
transition it is of the order of 1 cm. This is also the order of magnitude of the 
distance between the final quark-gluon plasma regions. Unfortunately, this is two 
orders of magnitude too small to affect nucleosynthesis. This result is rather robust 
against reasonable variations in any of the input parameters. 



6 Conclusion 

At temperatures exceeding approximately 160-180 MeV matter takes the form of 
quark-gluon plasma rather than hadronic gas. In principle QCD can be used to 
calculate the partition function and all other properties of interest. Perturbation 
theory has been employed to calculate the equation of state at asymptotically high 
temperatures to order g 5 ~ ct'J 2 , but the results diverge at low temperature and 
cannot predict the properties of a phase transition. For this purpose numerical 
Monte Carlo simulations of lattice gauge theory are necessary. For the real world 
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with 2 light flavors and one moderately light flavor of quark there is a big increase of 
the energy density in a narrow window of temperature centered on 170 MeV. Despite 
twenty years of effort it is still not known whether there is a true thermodynamic 
phase transition of first or second order or just a very rapid crossover. This is due 
to the great difficulty of computation with light dynamical quarks on the lattice. 

Nucleosynthesis is affected by remnant inhomogeneities in the baryon to entropy 
ratio and isospin if the high baryon density regions immediately following a QCD 
phase transition are separated by at least 1 m. A set of dynamical equations can be 
written and solved for the evolution of the universe through such a phase transition 
all the way to completion. The evolution of the temperature and hadronic volume 
fraction as functions of time and scale factor are hardly any different than the results 
of an idealized Maxwell construction. The information not available in the latter 
construction is the length scale of the inhomogeneities, that is, bubble sizes and 
so on. The characteristic distance between the last regions of quark-gluon plasma 
seem to be on the order of 1 cm, too small to affect nucleosynthesis. However, there 
are qualifications and improvements can be made. For example, when the fraction 
of space occupied by bubbles exceeds about 50%, interactions among the bubbles 
probably cannot be neglected. It is unlikely, though, that further improvements 
in the dynamics would qualitatively change the current picture of the transition. 
Indeed, crude estimates of the effects of bubble fusion on the dynamics of the QCD 
transition in heavy ion collisions indicate that the transition completes only a little 
faster, and that the average bubble size is greater p2 |. At least that is in the right 
direction to be interesting. 

The Relativistic Heavy Ion Collider (RHIC) is an accelerator at Brookhaven 
National Laboratory designed and constructed explicitly to search for quark-gluon 
plasma and the nature of a possible phase transition. Operations began in the sum- 
mer of 2000, and a wealth of data is expected soon. Essentially RHIC is attempting 
to reproduce the first microsecond of the early universe, to produce matter with 
temperatures exceeding well over 200 MeV. Many of us look forward to exploitation 



of this little big bang in the laboratory to cosmology [33 
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